Secreted novel AID/APOBEC-like deaminase 1 (SNAD1) - a new important player in fish immunology.

The AID/APOBECs are a group of zinc-dependent cytidine deaminases that catalyse the deamination of bases in nucleic acids, resulting in a cytidine to uridine transition. Secreted novel AID/APOBEC-like deaminases (SNADs), characterized by the presence of a signal peptide are unique among all of intracellular classical AID/APOBECs, which are the central part of antibody diversity and antiviral defense. To date, there is no available knowledge on SNADs including protein characterization, biochemical characteristics and catalytic activity. We used various in silico approaches to define the phylogeny of SNADs, their common structural features, and their potential structural variations in fish species. Our analysis provides strong evidence of the universal presence of SNAD1 proteins/transcripts in fish, in which expression commences after hatching and is highest in anatomical organs linked to the immune system. Moreover, we searched published fish data and identified previously, "uncharacterized proteins" and transcripts as SNAD1 sequences. Our review into immunological research suggests SNAD1 role in immune response to infection or immunization, and interactions with the intestinal microbiota. We also noted SNAD1 association with temperature acclimation, environmental pollution and sex-based expression differences, with females showing higher level. To validate in silico predictions we performed expression studies of several SNAD1 gene variants in carp, which revealed distinct patterns of responses under different conditions. Dual sensitivity to environmental and pathogenic stress highlights its importance in the fish and potentially enhancing thermotolerance and immune defense. Revealing the biological roles of SNADs represents an exciting new area of research related to the role of DNA and/or RNA editing in fish biology.

APOBEC2 safeguards skeletal muscle cell fate through binding chromatin and regulating transcription of non-muscle genes during myoblast differentiation.

The apolipoprotein B messenger RNA editing enzyme, catalytic polypeptide (APOBEC) family is composed of nucleic acid editors with roles ranging from antibody diversification to RNA editing. APOBEC2, a member of this family with an evolutionarily conserved nucleic acid-binding cytidine deaminase domain, has neither an established substrate nor function. Using a cellular model of muscle differentiation where APOBEC2 is inducibly expressed, we confirmed that APOBEC2 does not have the attributed molecular functions of the APOBEC family, such as RNA editing, DNA demethylation, and DNA mutation. Instead, we found that during muscle differentiation APOBEC2 occupied a specific motif within promoter regions; its removal from those regions resulted in transcriptional changes. Mechanistically, these changes reflect the direct interaction of APOBEC2 with histone deacetylase (HDAC) transcriptional corepressor complexes. We also found that APOBEC2 could bind DNA directly, in a sequence-specific fashion, suggesting that it functions as a recruiter of HDAC to specific genes whose promoters it occupies. These genes are normally suppressed during muscle cell differentiation, and their suppression may contribute to the safeguarding of muscle cell fate. Altogether, our results reveal a unique role for APOBEC2 within the APOBEC family.

Targeting APOBECs in cancer: It's about timing.

APOBEC3 cytidine deaminases have emerged as key drivers of mutagenesis in a wide spectrum of tumor types and are now appreciated to play a causal role in driving tumor evolution and drug resistance. As efforts to develop APOBEC3 inhibitors progress, understanding the timing and consequences of APOBEC3-mediated mutagenesis in distinct clinical contexts will be critical for guiding the development of anti-cancer therapeutic strategies.

Mesoscale DNA features impact APOBEC3A and APOBEC3B deaminase activity and shape tumor mutational landscapes.

Antiviral DNA cytosine deaminases APOBEC3A and APOBEC3B are major sources of mutations in cancer by catalyzing cytosine-to-uracil deamination. APOBEC3A preferentially targets single-stranded DNAs, with a noted affinity for DNA regions that adopt stem-loop secondary structures. However, the detailed substrate preferences of APOBEC3A and APOBEC3B have not been fully established, and the specific influence of the DNA sequence on APOBEC3A and APOBEC3B deaminase activity remains to be investigated. Here, we find that APOBEC3B also selectively targets DNA stem-loop structures, and they are distinct from those subjected to deamination by APOBEC3A. We develop Oligo-seq, an in vitro sequencing-based method to identify specific sequence contexts promoting APOBEC3A and APOBEC3B activity. Through this approach, we demonstrate that APOBEC3A and APOBEC3B deaminase activity is strongly regulated by specific sequences surrounding the targeted cytosine. Moreover, we identify the structural features of APOBEC3B and APOBEC3A responsible for their substrate preferences. Importantly, we determine that APOBEC3B-induced mutations in hairpin-forming sequences within tumor genomes differ from the DNA stem-loop sequences mutated by APOBEC3A. Together, our study provides evidence that APOBEC3A and APOBEC3B can generate distinct mutation landscapes in cancer genomes, driven by their unique substrate selectivity.

Distinguishing preferences of human APOBEC3A and APOBEC3B for cytosines in hairpin loops, and reflection of these preferences in APOBEC-signature cancer genome mutations.

The APOBEC3 enzymes convert cytosines in single-stranded DNA to uracils to protect against viruses and retrotransposons but can contribute to mutations that diversify tumors. To understand the mechanism of mutagenesis, we map the uracils resulting from expression of APOBEC3B or its catalytic carboxy-terminal domain (CTD) in Escherichia coli. Like APOBEC3A, the uracilomes of A3B and A3B-CTD show a preference to deaminate cytosines near transcription start sites and the lagging-strand replication templates and in hairpin loops. Both biochemical activities of the enzymes and genomic uracil distribution show that A3A prefers 3 nt loops the best, while A3B prefers 4 nt loops. Reanalysis of hairpin loop mutations in human tumors finds intrinsic characteristics of both the enzymes, with a much stronger contribution from A3A. We apply Hairpin Signatures 1 and 2, which define A3A and A3B preferences respectively and are orthogonal to published methods, to evaluate their contribution to human tumor mutations.

Engineered CBEs based on Macaca fascicularis A3A with improved properties for precise genome editing.

Cytidine deaminase defines the properties of cytosine base editors (CBEs) for C-to-T conversion. Replacing the cytidine deaminase rat APOBEC1 (rA1) in CBEs with a human APOBEC3A (hA3A) improves CBE properties. However, the potential CBE application of macaque A3A orthologs remains undetermined. Our current study develops and evaluates engineered CBEs based on Macaca fascicularis A3A (mA3A). Here, we demonstrate that BE4-mA3A and its RNA-editing-derived variants exhibit improved CBE properties, except for DNA off-target activity, compared to BE3-rA1 and BE4-rA1. Unexpectedly, deleting Ser-Val-Arg (SVR) in BE4-mA3A dramatically reduces DNA and RNA off-target activities and improves editing accuracy, with on-target efficiency unaffected. In contrast, a chimeric BE4-hA3A-SVR+ shows editing efficiency increased by about 50%, with other properties unaffected. Our findings demonstrate that mA3A-based CBEs could provide prototype options with advantages over rA1- and hA3A-based CBEs for further optimization, highlighting the importance of the SVR motif in defining CBE intrinsic properties.

Whole-Genome Mapping of Epigenetic Modification of 5-Formylcytosine at Single-Base Resolution by Chemical Labeling Enrichment and Deamination Sequencing.

DNA cytosine methylation (5-methylcytosine, 5mC) is a predominant epigenetic modification that plays a critical role in a variety of biological and pathological processes in mammals. In active DNA demethylation, the 10-11 translocation (TET) dioxygenases can sequentially oxidize 5mC to generate three modified forms of cytosine, 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). Beyond being a demethylation intermediate, recent studies have shown that 5fC has regulatory functions in gene expression and chromatin organization. While some methods have been developed to detect 5fC, genome-wide mapping of 5fC at base resolution is still highly desirable. Herein, we propose a chemical labeling enrichment and deamination sequencing (CLED-seq) method for detecting 5fC in genomic DNA at single-base resolution. The CLED-seq method utilizes selective labeling and enrichment of 5fC-containing DNA fragments, followed by deamination mediated by apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (APOBEC3A or A3A) and sequencing. In the CLED-seq process, while all C, 5mC, and 5hmC are interpreted as T during sequencing, 5fC is still read as C, enabling the precise detection of 5fC in DNA. Using the proposed CLED-seq method, we accomplished genome-wide mapping of 5fC in mouse embryonic stem cells. The mapping study revealed that promoter regions enriched with 5fC overlapped with H3K4me1, H3K4me3, and H3K27ac marks. These findings suggest a correlation between 5fC marks and active gene expression in mESCs. In conclusion, CLED-seq is a straightforward, bisulfite-free method that offers a valuable tool for detecting 5fC in genomes at a single-base resolution.

Activation-induced cytidine deaminase causes recurrent splicing mutations in diffuse large B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma. A major mutagenic process in DLBCL is aberrant somatic hypermutation (aSHM) by activation-induced cytidine deaminase (AID), which occurs preferentially at RCH/TW sequence motifs proximal to transcription start sites. Splice sequences are highly conserved, rich in RCH/TW motifs, and recurrently mutated in DLBCL. Therefore, we hypothesized that aSHM may cause recurrent splicing mutations in DLBCL. In a meta-cohort of > 1,800 DLBCLs, we found that 77.5% of splicing mutations in 29 recurrently mutated genes followed aSHM patterns. In addition, in whole-genome sequencing (WGS) data from 153 DLBCLs, proximal mutations in splice sequences, especially in donors, were significantly enriched in RCH/TW motifs (p < 0.01). We validated this enrichment in two additional DLBCL cohorts (N > 2,000; p < 0.0001) and confirmed its absence in 12 cancer types without aSHM (N > 6,300). Comparing sequencing data from mouse models with and without AID activity showed that the splice donor sequences were the top genomic feature enriched in AID-induced mutations (p < 0.0001). Finally, we observed that most AID-related splice site mutations are clonal within a sample, indicating that aSHM may cause early loss-of-function events in lymphomagenesis. Overall, these findings support that AID causes an overrepresentation of clonal splicing mutations in DLBCL.

Monitoring APOBEC3A protein levels in human cancer cells.

The APOBEC3 family of cytosine deaminases, which target single-stranded DNA and RNA of viruses and retroelements as part of the innate immune defense, generate mutations in many human cancers. Although the APOBEC3A paralog is a major endogenous source of these mutations, low APOBEC3A mRNA levels and protein abundance have hampered functional characterization. Extensive homology across APOBEC3 paralogs have further challenged the development of specific detection reagents. Here, we describe the isolation and use of monoclonal antibodies with specificity for APOBEC3A and the APOBEC3A/APOBEC3B/APOBEC3G proteins. We provide protocols and technical advice for detection and measurement of APOBEC3A protein across human cancer cell lines using standard immunoblotting and immunofluorescence protocols.

Compact zinc finger architecture utilizing toxin-derived cytidine deaminases for highly efficient base editing in human cells.

Nucleobase editors represent an emerging technology that enables precise single-base edits to the genomes of eukaryotic cells. Most nucleobase editors use deaminase domains that act upon single-stranded DNA and require RNA-guided proteins such as Cas9 to unwind the DNA prior to editing. However, the most recent class of base editors utilizes a deaminase domain, DddAtox, that can act upon double-stranded DNA. Here, we target DddAtox fragments and a FokI-based nickase to the human CIITA gene by fusing these domains to arrays of engineered zinc fingers (ZFs). We also identify a broad variety of Toxin-Derived Deaminases (TDDs) orthologous to DddAtox that allow us to fine-tune properties such as targeting density and specificity. TDD-derived ZF base editors enable up to 73% base editing in T cells with good cell viability and favorable specificity.

Identification of significant m6A regulators and immune microenvironment characterization in ischemic stroke.

The role of m6A modification in the regulation of the immune microenvironment (IME) of ischemic stroke (IS) is barely known. Thus, we aim to investigate the impact of m6A modification on the IME of IS and its diagnostic value in IS. We comprehensively assessed the m6A modification patterns, the relationship between these modification patterns and the characteristics of the IME. The m6A modification patterns of individual IS sample were quantified by m6Ascore. The performance of m6A phenotype-related genes as potential biomarkers was evaluated by the area under the receiver operating characteristic curve. Experimental validation was also performed by qRT-PCR. Six dysregulated m6A regulators were identified and a classification model consisting of four key m6A regulators (METLL3, RBMX, RBM15B, YTDHF3) could distinguish IS and healthy control samples well. METTL3 and YTHDF3 are closely related to circulating neutrophil abundance. Two distinct m6A modification patterns were determined which differed in immunocyte abundance. We also identified six m6A phenotype-related genes (APOBEC3A, PTMA, FCGR3A, LOC440926, LOC649946, and FTH1L11), and further explored their biological function. Among them, APOBEC3A, FCGR3A, and FTH1L11 were positively associated with neutrophil abundance. APOBEC3A and FCGR3A were stable diagnostic m6A-associated genes in both the discovery and validation cohorts. This study reveals that m6A modification plays a non-negligible role in the formation of a diversified and complex IME in IS. The m6A phenotype-related genes could be diagnostic biomarkers of IS.

A Novel Heterozygous Variant in AICDA Impairs Ig Class Switching and Somatic Hypermutation in Human B Cells and is Associated with Autosomal Dominant HIGM2 Syndrome.

B cells and their secreted antibodies are fundamental for host-defense against pathogens. The generation of high-affinity class switched antibodies results from both somatic hypermutation (SHM) of the immunoglobulin (Ig) variable region genes of the B-cell receptor and class switch recombination (CSR) which alters the Ig heavy chain constant region. Both of these processes are initiated by the enzyme activation-induced cytidine deaminase (AID), encoded by AICDA. Deleterious variants in AICDA are causal of hyper-IgM syndrome type 2 (HIGM2), a B-cell intrinsic primary immunodeficiency characterised by recurrent infections and low serum IgG and IgA levels. Biallelic variants affecting exons 2, 3 or 4 of AICDA have been identified that impair both CSR and SHM in patients with autosomal recessive HIGM2. Interestingly, B cells from patients with autosomal dominant HIGM2, caused by heterozygous variants (V186X, R190X) located in AICDA exon 5 encoding the nuclear export signal (NES) domain, show abolished CSR but variable SHM. We herein report the immunological and functional phenotype of two related patients presenting with common variable immunodeficiency who were found to have a novel heterozygous variant in AICDA (L189X). This variant led to a truncated AID protein lacking the last 10 amino acids of the NES at the C-terminal domain. Interestingly, patients' B cells carrying the L189X variant exhibited not only greatly impaired CSR but also SHM in vivo, as well as CSR and production of IgG and IgA in vitro. Our findings demonstrate that the NES domain of AID can be essential for SHM, as well as for CSR, thereby refining the correlation between AICDA genotype and SHM phenotype as well as broadening our understanding of the pathophysiology of HIGM disorders.

Mayaro virus infection elicits a robust pro-inflammatory and antiviral response in human macrophages.

Mayaro virus (MAYV), the etiological agent of Mayaro fever (MAYF), is an emergent arbovirus pathogen belonging to Togaviridae family. MAYF is characterized by high inflammatory component that can cause long-lasting arthralgia that persists for months. Macrophages are viral targets and reservoirs, key components of innate immunity and host response. Given the importance of this pathogen, our aim was to determine the inflammatory and antiviral response of human monocyte-derived macrophages (MDMs) infected with MAYV. First, we established the replication kinetics of the virus. Thereafter, we determined the expression of pattern recognition receptors, NF-ĸB complex, interferons (IFNs), two interleukin 27 (IL27) subunits, IFN-stimulated genes (ISGs), and the production of cytokines/chemokines. We found that human MDMs are susceptible to MAYV infection in vitro, with a peak of viral particles released between 24- and 48-hours post-infection (h.p.i) at MOI 0.5, and between 12 and 24 h.p.i at MOI 1. Interestingly, we observed a significant decline in the production of infectious viral particles at 72 h.p.i that was associated with the induction of antiviral response and high cytotoxic effect of MAYV infection in MDMs. We observed modulation of several genes after MAYV infection, as well, we noted the activation of antiviral detection and response pathways (Toll-like receptors, RIG-I/MDA5, and PKR) at 48 h.p.i but not at 6 h.p.i. Furthermore, MAYV-infected macrophages express high levels of the three types of IFNs and the two IL27 subunits at 48 h.p.i. Moreover, we found higher production of IL6, IL1ß, CXCL8/IL8, CCL2, and CCL5 at 48 h.p.i as compared to 6 h.p.i. A robust antiviral response (ISG15, APOBEC3A, IFITM1, and MX2) was observed at 48 but not at 6 h.p.i. The innate and antiviral responses of MAYV-infected MDMs differ at 6 and 48 h.p.i. We conclude that MAYV infection induces robust pro-inflammatory and antiviral responses in human primary macrophages.

Discovery of cytosine deaminases enables base-resolution methylome mapping using a single enzyme.

Deaminases have important uses in modification detection and genome editing. However, the range of applications is limited by the small number of characterized enzymes. To expand the toolkit of deaminases, we developed an in vitro approach that bypasses a major hurdle with their toxicity in cells. We assayed 175 putative cytosine deaminases on a variety of substrates and found a broad range of activity on double- and single-stranded DNA in various sequence contexts, including CpG-specific deaminases and enzymes without sequence preference. We also characterized enzyme selectivity across six DNA modifications and reported enzymes that do not deaminate modified cytosines. The detailed analysis of diverse deaminases opens new avenues for biotechnological and medical applications. As a demonstration, we developed SEM-seq, a non-destructive single-enzyme methylation sequencing method using a modification-sensitive double-stranded DNA deaminase. The streamlined protocol enables accurate, base-resolution methylome mapping of scarce biological material, including cell-free DNA and 10 pg input DNA.

Stage II oesophageal carcinoma: peril in disguise associated with cellular reprogramming and oncogenesis regulated by pseudogenes.

INTRODUCTION: Pseudogenes have been implicated for their role in regulating cellular differentiation and organismal development. However, their role in promoting cancer-associated differentiation has not been well-studied. This study explores the tumour landscape of oesophageal carcinoma to identify pseudogenes that may regulate events of differentiation to promote oncogenic transformation. MATERIALS AND METHOD: De-regulated differentiation-associated pseudogenes were identified using DeSeq2 followed by 'InteractiVenn' analysis to identify their expression pattern. Gene expression dependent and independent enrichment analyses were performed with GSEA and ShinyGO, respectively, followed by quantification of cellular reprogramming, extent of differentiation and pleiotropy using three unique metrics. Stage-specific gene regulatory networks using Bayesian Network Splitting Average were generated, followed by network topology analysis. MEME, STREME and Tomtom were employed to identify transcription factors and miRNAs that play a regulatory role downstream of pseudogenes to initiate cellular reprogramming and further promote oncogenic transformation. The patient samples were stratified based on the expression pattern of pseudogenes, followed by GSEA, mutation analysis and survival analysis using GSEA, MAF and 'survminer', respectively. RESULTS: Pseudogenes display a unique stage-wise expression pattern that characterizes stage II (SII) ESCA with a high rate of cellular reprogramming, degree of differentiation and pleiotropy. Gene regulatory network and associated topology indicate high robustness, thus validating high pleiotropy observed for SII. Pseudogene-regulated expression of SOX2, FEV, PRRX1 and TFAP2A in SII may modulate cellular reprogramming and promote oncogenesis. Additionally, patient stratification-based mutational analysis in SII signifies APOBEC3A (A3A) as a potential hallmark of homeostatic mutational events of reprogrammed cells which in addition to de-regulated APOBEC3G leads to distinct events of hypermutations. Further enrichment analysis for both cohorts revealed the critical role of combinatorial expression of pseudogenes in cellular reprogramming. Finally, survival analysis reveals distinct genes that promote poor prognosis in SII ESCA and patient-stratified cohorts, thus providing valuable prognostic bio-markers along with markers of differentiation and oncogenesis for distinct landscapes of pseudogene expression. CONCLUSION: Pseudogenes associated with the events of differentiation potentially aid in the initiation of cellular reprogramming to facilitate oncogenic transformation, especially during SII ESCA. Despite a better overall survival of SII, patient stratification reveals combinatorial de-regulation of pseudogenes as a notable marker for a high degree of cellular differentiation with a unique mutational landscape.

DNA flexibility can shape the preferential hypermutation of antibody genes.

Antibody-coding genes accumulate somatic mutations to achieve antibody affinity maturation. Genetic dissection using various mouse models has shown that intrinsic hypermutations occur preferentially and are predisposed in the DNA region encoding antigen-contacting residues. The molecular basis of nonrandom/preferential mutations is a long-sought question in the field. Here, we summarize recent findings on how single-strand (ss)DNA flexibility facilitates activation-induced cytidine deaminase (AID) activity and fine-tunes the mutation rates at a mesoscale within the antibody variable domain exon. We propose that antibody coding sequences are selected based on mutability during the evolution of adaptive immunity and that DNA mechanics play a noncoding role in the genome. The mechanics code may also determine other cellular DNA metabolism processes, which awaits future investigation.

Cooperativity between Cas9 and hyperactive AID establishes broad and diversifying mutational footprints in base editors.

The partnership of DNA deaminase enzymes with CRISPR-Cas nucleases is now a well-established method to enable targeted genomic base editing. However, an understanding of how Cas9 and DNA deaminases collaborate to shape base editor (BE) outcomes has been lacking. Here, we support a novel mechanistic model of base editing by deriving a range of hyperactive activation-induced deaminase (AID) base editors (hBEs) and exploiting their characteristic diversifying activity. Our model involves multiple layers of previously underappreciated cooperativity in BE steps including: (i) Cas9 binding can potentially expose both DNA strands for 'capture' by the deaminase, a feature that is enhanced by guide RNA mismatches; (ii) after strand capture, the intrinsic activity of the DNA deaminase can tune window size and base editing efficiency; (iii) Cas9 defines the boundaries of editing on each strand, with deamination blocked by Cas9 binding to either the PAM or the protospacer and (iv) non-canonical edits on the guide RNA bound strand can be further elicited by changing which strand is nicked by Cas9. Leveraging insights from our mechanistic model, we create novel hBEs that can remarkably generate simultaneous C > T and G > A transitions over >65 bp with significant potential for targeted gene diversification.

Cytidine Deaminase Resolves Replicative Stress and Protects Pancreatic Cancer from DNA-Targeting Drugs.

Cytidine deaminase (CDA) functions in the pyrimidine salvage pathway for DNA and RNA syntheses and has been shown to protect cancer cells from deoxycytidine-based chemotherapies. In this study, we observed that CDA was overexpressed in pancreatic adenocarcinoma from patients at baseline and was essential for experimental tumor growth. Mechanistic investigations revealed that CDA localized to replication forks where it increased replication speed, improved replication fork restart efficiency, reduced endogenous replication stress, minimized DNA breaks, and regulated genetic stability during DNA replication. In cellular pancreatic cancer models, high CDA expression correlated with resistance to DNA-damaging agents. Silencing CDA in patient-derived primary cultures in vitro and in orthotopic xenografts in vivo increased replication stress and sensitized pancreatic adenocarcinoma cells to oxaliplatin. This study sheds light on the role of CDA in pancreatic adenocarcinoma, offering insights into how this tumor type modulates replication stress. These findings suggest that CDA expression could potentially predict therapeutic efficacy and that targeting CDA induces intolerable levels of replication stress in cancer cells, particularly when combined with DNA-targeted therapies. SIGNIFICANCE: Cytidine deaminase reduces replication stress and regulates DNA replication to confer resistance to DNA-damaging drugs in pancreatic cancer, unveiling a molecular vulnerability that could enhance treatment response.

A computational exploration of global and temporal dynamics of selection pressure on HIV-1 Vif polymorphism.

Virion infectivity factor (Vif), an accessory protein of HIV-1 (human immunodeficiency virus type 1), antagonizes host APOBEC3 protein (apolipoprotein B mRNA editing enzyme, catalytic polypeptide 3) or A3 via proteasomal degradation, facilitating viral replication. HLA (Human leukocyte antigens) alleles, host restriction factors, and error-prone reverse transcription contribute to the global polymorphic dynamics of HIV, impacting effective vaccine design. Our computational analysis of over 50,000 HIV-1 M vif sequences from the Los Alamos National Laboratory (LANL) database (1998-2021) revealed positive selection pressure on the vif gene (nonsynonymous to synonymous ratio, dn/ds=1.58) and an average entropy score of 0.372 in protein level. Interestingly, over the years (1998-2021), a decreasing trend of dn/ds (1.68 to 1.47) and an increasing trend of entropy (0.309 to 0.399) was observed. The predicted mutational frequency against Vif consensus sequence decreased over time (slope = -0.00024, p < 0.0001). Sequence conservation was observed in Vif functional motifs F1, F2, F3, G, BC box, and CBF ß binding region, while variability was observed mainly in N- and C- terminal and Zinc finger region, which were dominantly under immune pressure by host HLA-I-restricted CD8+ T cell. Computational analysis of ∆∆Gstability through protein stability prediction tools suggested that missense mutation may affect Vif stability, especially in the Vif-A3 binding interface. Notably, mutations R17K and Y44F in F1 and G box were predicted to destabilize the Vif-A3 binding interface by altering bond formations with adjacent amino acids. Therefore, our analysis demonstrates Vif adaptation with host physiology by maintaining sequence conservation, especially in A3 interacting functional motifs, highlighting important therapeutic candidate regions of Vif against HIV-1 infections.